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Abstract 

Recent observational indications of an accelerating universe enhance the interest in 
studying models with a cosmological constant. We investigate cosmological expan- 
sion (FRW metric) with A > for a general linear equation of state p — wp, w > — 1, 
so that the interplay between cosmological vacuum and quintessence is allowed, as 
well. 

Four closed-form solutions (flat universe with any w, and w = 1/3, —1/3, —2/3) 
are given, in a proper compact representation. Various estimates of the expansion 
are presented in a general case when no closed-form solutions are available. For 
the open universe a simple relation between solutions with different parameters is 
established: it turns out that a solution with some w and (properly scaled) A is 
expressed algebraically via another solution with special different values of these 
parameters. 

The expansion becomes exponential at large times, and the amplitude at the 
exponent depends on the parameters. We study this dependence in detail, deriving 
various representations for the amplitude in terms of integrals and series. The closed- 
form solutions serve as benchmarks, and the solution transformation property noted 
above serves as a useful tool. Among the results obtained, one is that for the open 
universe with relatively small cosmological constant the amplitude is independent 
of the equation of state. Also, estimates of the cosmic age through the observable 
ratio Qji/Qjrf and parameter w are given; when inverted, they provide an estimate 
of w, i. e., the state equation, through the known and age of the universe. 
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1 Introduction 



Recent data on the brightness of distant SN la ([1,2]; see [3] and the references therein), 
as well as evidence coming from the cosmic age, large scale structure, and cosmic mi- 
crowave background anisotropy combined with the cluster dynamics, indicate, most 
probably, that the observed cosmological expansion is accelerating. Perhaps the most 
natural, although definitely not unique, reason for this acceleration is the presence of a 
cosmic vacuum of nonzero energy and pressure; because of that, investigation of vari- 
ous cosmological models including a positive cosmological constant becomes interesting, 
once again. In this paper we study one such model involving matter with an arbitrary 
(linear) equation of state; in particular, an interplay between vacuum and quintessence 
[4] is studied. 

We consider a Friedmann cosmology described by the Friedmann-Robertson- Walker 
metric , 

ds 2 = _ dt 2 + 2 (t) J^_A!_ + r 2 do 2^ ^ (1) 

where dQ is the element of solid angle and k = —1,0, or 1, according to whether the 
universe is open, flat, or closed (we use the units with G = c = 1). As usual, the energy- 
momentum tensor, consistent with homogeneity and isotropy of the universe, is the one 
corresponding to a perfect fluid, described by the energy density p and pressure p. The 
latter are assumed to satisfy the linear equation of state 

p = wp, w = const > — 1 . (2) 

The range of values of w is chosen on the grounds that w = —1 corresponds to the 
vacuum equation of state, and vacuum contributions are already included in our model 
with the cosmological constant A > 0. Moreover, no physical reasons seem to be known 
so far to justify the values of w which are more negative that the vacuum one, w = — 1, 
and the cosmological evolution in this range is very different. On the other hadn, it is 
not easy to imagine the kind of physics that would lead to matter with the pressure 
larger than its density, but the following analysis does not depend at all on whether w 
is larger or smaller than one, so we do not limit this parameter from above. Condition 
(|2|) allows for quintessence , that is, for — 1 < w < 0; investigation of its influence on 
cosmological expansion is one of the goals of this paper. 
Under these assumptions, Einstein's equations reduce to 

lda \ 2 o i a q k dp o P+P da ro\ 
= 8irp + A - 3 — , — = -3 — , (3 

a at J a z at a at 

the second one expressing energy conservation. By (|2|), this equation is easily integrated 
giving the density in terms of the scale factor, 

M 

P O 3(to+1) ' W 

where the constant M > characterizes the abundance of matter in the universe. A 
substituion of (||) in the first of equations @ produces the equation for the scale factor 
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da\ 2 8vrM 1 A » , 

+ — a — k ; (5) 



a(t) only, 

{.dtj 3 a 3w+1 ' 3 

we append it with the Big Bang initial condition 

o(0) = . (6) 

As seen from (|j), the initial value of the density, p(0), is then infinite for w > —1/3, and 
zero for — 1 < w < —1/3; it is finite and positive for w = —1/3, which value corresponds 
to the so called Einstein quintessence^]. 

Equation (|5|) with the initial condition (^) constitute the problem studied below; its 
solution, together with formula (|j), specifies completely the corresponding cosmological 
expansion at all times. 

2 Large— Time Expression for the Scale Factor 

2.1 Problem Reformulation and Parameter Range 

For convenience, we rescale the variables according to 

[3 , . /8ttM\ 1 / 3 ^ +1 ) . . 

* = Va t ' a ^ = {-jr) y(r) ' (7) 



and introduce new parameters: 

q / A \ 2/3(to + l) 

v = Bw + l, u>-2; u = 4 f — ) = 3 (87rM)- 2 ^+ 2 ) ( A )-^+ 2 ) > . 

(8) 

Denoting by a dot the derivative with respect to r, we rewrite (||) and (||) as 

y 2 = y~ u + y 2 -ha, y(0) = 0; i/ > -2; fc = 0,±l; w>0; (9) 

since we are studying the expansion, a positive square root of the right hand side of (||) 
is, of course, always taken for y. 

The equation and initial condition (^) are consistent for any value of parameter kui 
if v is positive, because at small times y 2 ~ y~ v > 0. For v = at small times one gets 
y 2 ~ 1 — ku>, so u) is limited by unity for the closed universe (k = 1). When —2 < v < 0, 
the governing equation in the beginning of the expansion becomes y 2 ~ — ku>, so only 
open and flat universe (k = — 1, 0) solutions are possible. 

Moreover, we are interested only in the solutions describing an unbounded cosmo- 
logical expansion, i. e., those with y(r) — » oo at r — > cxd. A simple mechanical analogy 
helps to establish parameter limits ensuring that. Namely, equation (S) describes the 
motion of a particle of a unit mass along the positive y axis, starting at the origin, in 
the potential 

V(y) = - ( y~ u + y 2 ) < ; 



1 The term 'Einstein quintessence' was coined recently in [5] on the grounds of the effective equation 
of state for Einstein's cosmological solution [6] of 1917. 
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the total energy of the particle is — ku. When —2 < v < 0, the potential is monotonically 
decreasing, with the largest value V(0) = 0. Hence any such motion with a non-negative 
total energy (k = —1,0, open and flat universe) is infinite. The same monotonicity of 
the potential holds for v = 0, but the largest value now is V(0) = —1, which yields 
the already known restriction u < 1 for any infinite motion in the case k = 1 (closed 
universe) . 

Finally, for v > the potential is strictly negative and non-monotonic: at y = y v > 
such that V (y v ) = it reaches a unique maximum V(y v ) = —uj u given by 

^ = min (y- v + y 2 ) = y^* + jfi = y v = (0 ^ + } . (10) 

An infinite motion requires the total energy to be larger than that maximum, —kui> 
—u u , which is automatically fulfilled for the open and flat universe; however, for the 
closed universe (k = 1), the parameter u> proves to be limited from above, to < ui u . This 
constraint has a clear physical meaning: to overwhelm the curvature effect in the closed 
universe and get an unbounded expansion, one has to have enough vacuum repulsion, 
that is, a large enough A. But for v > the parameter to is proportional to a negative 
power of A [see @], hence u> cannot be too large. 

Note that y u — > +0 and oj u — > 1 + when v — > +0, which allows us to describe the 
cases of positive and zero v uniformly. Also, both y u and ui u tend to unity from above 
when v — > oo, hence they have at least one maximum. In fact, each of the functions, 
plotted in Fig.l, has one maximum exactly: maxuv = 2 at v = 2 (for some reason, 
it is relativistic gas that has the largest 'store of expansion' in the closed universe), 
and maxy u ~ 1.149 at v ~ 7.182. Some values of particular interest are: for v = 
{w = —1/3, Einstein quintessence), y u = 0, uj u = 1; for u = 1 (w = 0, pressureless 
'dust'), y v = 2- 1 / 3 ps 0.794, lo u = 3 x 2~ 2 / 3 w 1.890; for v = 2 (w = 1/3, relativistic 
gas), y v = 1, oj v = 2; for v = 4 (w = 1, Zeldovich's superstiff fluid), y v = 2 1 / 6 ss 1.122, 
and again w„ = 3x 2~ 2 / 3 « 1.890. 

Summarizing, the range of parameters of interest is 

k = 1 : v > 0, < to < uj v \ A; = 0,-1: v > -2, < uj < oo . (11) 
2.2 Expansion Amplitude and Its Basic Properties 

With parameters in the range (0), the cosmological expansion goes indefinitely [unique 
solution to problem (^) grows unboundedly], and at large times it is exponential in time, 
since for y S> max{l, uj} the equation becomes y ~ y [for k = 0, 1 the large time 
condition is apparently relaxed to just y S> 1]. Therefore 

y(r) = y(r, u, kuj) = A u (kuj) e T [1 + o(l)] , r -» oo , (12) 

and we are particularly interested in the behavior of the expansion amplitude A v {kui) 
for the whole parameter range (|ll]). For v > and kui = u) u (i. e., k = 1, u; = w^) the 
solution y(r) reaches the (unstable) rest point y v in an infinite time, so A u (oj u ) should 
be zero, which remains to be proved. 
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We derive an explicit expression for the amplitude by integrating equation (0), 
r(y) = r(y, v, ku) 1 



\fx v + x 2 — kuj 



(13) 



For y, t — ► oo the integral here diverges logarithmically, as it should. To get the desired 
expression for the amplitude, one needs to extract this logarithm of y, that is, to add 
and subtract something in the integrand allowing for a closed-form quadrature. The 
most natural choice seems the same integrand with k = 0, in which case one writes: 



r(y) 



In 



In 



u + 2 



2 2/(, + 2) y 



y{«+W + ^1 + y(^+2) 



+ 



1 



1 



+ 



Vx~ 
1 



+ x^ 



+ x 1 



v + x 2 — kuj \fx v + x 
dx + o(l), y 



dx 



(14) 



[unlike (|l3|), the integral here converges at infinity]. This expression can be combined 
with the equality ( |l~2| ) to write the desired fomula for the amplitude as 



A u (kLj) = A v (0)e-^ ku \ 
where 



1 



1 



\/^ u + x 2 — kuj va 



+ 



(15) 



A v (0) = 2- 2 '^ 



(note that, naturally, (f> u (0) = 0). 

Two features of the expansion amplitude behavior can be established immediately. 
First, when koj = u> — > oj v — 0, the integral (|l5|) for (j) v (ku) diverges at x = so <p v {koj) 
goes to +oo; therefore A v {uj u ) = 0, as expected. 

Next, we have, for a given v: 



din A u (kuj) 



d(fi u (ku>) 



d(kiv) d(kuj) 
and the latter derivative is readily found from (|l~5|): 



d(j) v {kuj) 
d{ku>) 



dx 



> 



(16) 



(17) 



lo {x~ v + x 2 - kuj) 3/2 

Therefore A u (kto) is a decreasing function of kuj for any given relevant v\ in particular, 

A u {-lo) > A„{0)[> A u (uj), v > 0, u < u v ] . (18) 

That means that, as one expects intuitively, for the same parameters w, M, and A, the 
amplitude for the open universe is larger than that for the flat universe which, in its 
turn, is larger than the amplitude for the closed universe if the latter exists. Generally, 
by figuring out the signs of the curvature and parameter v and taking into account 
expression (g) for u, one can see that, given k = dbl, u, and the matter abundance M, 
the growth amplitude is larger for larger cosmological constant A. That clearly makes 
physical sense: the larger the contribution of the vacuum repulsion, the stronger the 
expansion is. 
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3 Closed— Form Cosmological Solutions 



The closed-form and, especially, elementary solutions are always most valuable; in the 
context of this study, they specifically shed light on how the expansion amplitude be- 
haves. 

It is not difficult to see that all the integrals involved in the general cosmological 
solution of the previous section prove to be elementary in the following four cases. 



3.1 Flat Universe: k = 



For k = equation (12) gives the result presented in part in the previous section; here 
is how it looks in full: 



T (y) 



u + 2 



In 



/"+2)/2 + ^1 + y{u+2) 



y( T ) 



sinh 



u + 2 \-|2/("+2) 



• (19) 



Thus for flat universe the amplitude of the exponential growth of the scale factor is 

A v (0) = 2- 2 /<r+ 2 \ 0<A,(0)<1. (20) 
Note also the small time asymptotics expression, 

V + 2 \ 2 /(^+2) 



y{ T > 



[l + o(l)], 



which is valid for all the values of koo and v > (the curvature and cosmological constant 
do not play any role in the early universe). For v = 0, and v < 0, k = —1, the expansion 
becomes linear with the time, in the beginning. 



3.2 Relativistic Gas: w = 1/3 

Here v = 2, and we have: 



, , 1 , 2y 2 - kuj + 2^y 4 - kujy 2 + 1 
T{V) = 2 b 2—k^ 



' 2 — kuj , 2 + koj , . kuj 



y(r) = y — - — exp(2r) — exp(— 2r) + — = y sinh(2r) — kuj sinh r , (21) 

so that 

, /, n V2 — kuj 

A 2 {ku) = (22) 

(all these results agree with the above ones for k = 0). Note that UJ2 = 2 [the definition 
of uj u is given in (|l^) ] , and the expression (^2|) for the amplitude is well defined exactly 
up to uj = L02 = 2, for the closed universe (k = 1). 
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3.3 Einstein quintessence: w = — 1/3 

In this case we have v = 0, and 

r(y) = In - — , y(r) = \/l — kuj sinh r , (23) 



\/l — kuj 



hence 



Ao(ku;) = , (24) 



again in agreement with section 3.1. Of course, as seen from fllQ|), ojq = 1 

3.4 Quintessence with w = — 2/3 

Finally, here u = —1, thus 



r(,) = ln^ + 2 ^!_ + ^, . = 0,-1; 



y(r) = — — — — exp(r) + - [exp(— r) — 1] = \/—kuj sinh r + sinh 2 (t/2) , (25) 

and 

1 + 2y/-k0J . . 

A_i(A:w) = , (26) 

Naturally, these expressions remain physically meaningful only for = — 1, 0, coinciding 
with the results of section 3.1 in the latter case of flat universe. 

Note that all the above solutions are known (see c. f. [7,8]), although not always 
in exactly the same form. In the proper current theoretical context, the state equation 
with uu = —1/3 is called the equation of state for strings, while w = —2/3 provides the 
equation of state for domain walls (e. g. [9]). 

3.5 Closed— Form Solutions in Terms of Elliptic Functions 

By definition, an elliptic integral is the integral of any rational function of the integration 
variable and the square root of a polynomial of this variable whose degree is three or 
four. Any elliptic integral can be reduced to a combination of elementary functions and 
three standarad elliptic integrals (see c. f. [10]). 



The form of integrals (13) and ( [lq ) determining our cosmological solution and its 
large-time asymptotics prompts an arbitrary power change of the integration variable to 
try reducing them to elliptic integrals. Such an attempt turns out to be successful for 
four different positive values of parameter v. By the transformation property of section 
5, four reciprocal negative values should be added to the list, making a total of eight 
different exact solutions in terms of elliptic functions. Since these representations are 
rather cumbersome and not very informative, we choose not to give them here explicitly; 
instead, in Table 1 below we show the corresponding values of v and w, and the proper 
substitution making the integrals elliptic. 
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Table 1. Elliptic function solutions for v > 



V 


u> 


Substitution 


6 


5/3 


x = e /2 


4 


1 


x = r 1/2 


1 





X = f" 1 


2/3 


-1/9 





Note that u> = describes pressureless matter ('dust'), and w = 1 corresponds to 
Zeldovich's superstiff fluid [11]. According to (|39[), the reciprocal negative values of f 
allowing for the elliptic function solutions are v = —3/2 — 4/3, —2/3, —1/2, k = —1. 



4 Estimates of the Scale Factor and Expansion Amplitude 

In the general case, no closed-form solution is available, so it is useful to have at least 
some estimates illustrating the behavior of cosmological expansion. Here we derive two 
groups of such estimates differing by the parameters involved in them. It is not difficult 
to check that all the exact elementary solutions satisfy these estimates, as they must. 



4.1 Estimates Independent of the State Equation 

Some estimates of the scale factor and the expansion amplitude follow immediately from 
(M). Indeed, we have 



y=\Jy u + y 2 - kuj > \Jy 2 - kuj 
for all y > when k = —1, 0, and for y > y/co when k = 1. From here we obtain: 



— ^y 2 - kuj^j > 1, y > 0, k = -1; y > y/Z, k = l. (27) 

In the case of an open universe we integrate this between zero and an arbitrary moment 
of time, which, due to the initial condition (||), results in 



y + \jy 2 +lu > yfuj exp r . 

Solving this inequality for y provides a simple lower bound of the open universe expansion 
at all times: 

y(r) = y(r, u, — u) > y/u sinhr, r>0 (resp., y > 0), , k = — 1; . (28) 
Sending r — ► oo and comparing to (|i"2|), we find the estimate for the growth amplitude: 

A v (-u) > y/u/2, k = -1 . (29) 



In the case of closed universe we integrate inequality (|27| ) up to an arbitrary moment 
of time starting with [see (O)] 

Tu = T(y/u,v,u;) = / = , (30) 

Jo \J X v + x z — CO 
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with y changing, respectively, from y/U to y{j). Thus we arrive at 



y + \fy 2 -u > exp(r - r u ) ; 
solution of this inequality yields the estimate 

y (r) = y(r, v,uj) > cosh(r - t u ), t>t u (resp., y > y/u), k = 1 . (31) 



Again, in the limit r — > oo this inequality combined with ([12]) allows for the correspond- 
ing estimate of the expansion amplitude: 



A v (uj) > exp(-T w ), k = 1, v > . 



(32) 



Note that, by (fftj), t w — > oo when — > lo v — 0, hence the estimate (p2| ) does not 
contradict the amplitude approaching zero at u v . Clearly, all these estimates do not 
explicitly depend on v, except the dependence on it of r w in the case of a closed universe. 



4.2 Estimates Independent of Parameter uj 

On the physical grounds, one expects that the time in which the expansion reaches a 
given value of the scale factor, y, increases from open to flat to closed universe, 

r(y,u,-u) < r(y,i/,0)[< r(y,u,u), v > 0] , 

because of the curvature effect. Formally, this is easily derived from the general solution 
(fl3|) by dropping the term ku in the denominator of the integrand there, and choosing 
the inequality sign accordingly. Using the exact solution ( |T9| ) for the flat universe, we 
thus obtain: 



r < 



u + 2 



In 



j,("+2)/2 + ^/l + y (u+2) 



T > 



u + 2 



In 



y (H-2)/2 + .^1 + y (,+2) 



-1, ^> -2; 



1, i/ > 



(33) 



By inverting these, we get the estimates for the scale factor: 

2/(i/+2) 



y(r) > 



sinh ( — — — t 



k 



-1, ^> -2; 



(34) 



sinh ( — — — r 



2/(u+2) 



k = 1, i/ > . 



Naturally, the estimates for the growth amplitude implied by (|34j) are nothing else than 
(fl|), in view of @). 
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5 Transformation Property of Cosmological Solutions and 
Expansion Amplitude for the Open Universe 

5.1 Derivation of Transformation Property 

A family of cosmological solutions describing an open universe turns out to have a useful 
transformation property. It relates, in a rather simple algebraic way, a solution with a 
given set of parameters v, to to another with different values of these parameters. 

To establish this property, let us consider the initial value problem (|9|) with k = — 1, 
kuj = -co, and introduce a new unknown function z = z(r) = z(t,i/,—uj) instead of 
U = y( T ) = y{ T , u -> ~ UJ ) according to 



y = w l/2^/(H-2). 



2a; 1 / 2 
v + 2 



Since v + 2 > 0, the initial condition for z is the same as for y, 

2(0) = . 

On the other hand, equation @ after substitution (|35| ) becomes 



(35) 



(36) 



+ 

Introducing the new time, 



$■ = Co-/ 2 - 1 



+ z 



(4+2v)/(v+2) + z -(-2v/(v+2)) 



v + 2 

T = ~ T 



and denoting by a prime the derivative, we rewrite the last equation as 

, 2 



-(-2v/(u+2)) 



+ Z" +LO 



-0+2)/2 



(37) 



(38) 



Clearly, the initial value problem (38), (|3q ) coincides with the original problem (||) 
in which r is replaced with f, v with —2v/{y + 2), and u with that, by the 

uniqueness theorem, 

z = z(f, v, -u) =y((u + 2)r/2, -2i//(i/ + 2), -c^+ 2 )/ 2 

But then, according to (]35|) and (|3 

y(r, i/, -w) = a; 1 / 2 [y ((„ + 2)r/2, -2i//(i/ + 2), -^-^+ 2 '/ 2 ) 

(39) 

where the identity sign stresses that the equality holds for all moments of time, as well 
as for any u > — 2 and cj > 0. Sending r to infinity provides, by (12), the corresponding 
transformation property of the expansion amplitude, 



2/(z/+2) 



AJ-u) 
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1/2 



■2u/(u+2) [-V 



-(v+2)/2 



2/(u+2) 



k 



1 



(40) 
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In particular, our previous results (23) and fl2"4|) for v = satisfy (|39| ) and (§Q), respec- 
tively. 

Note also the corresponding transformation of cosmic ages as functions of the scale 
factor y, 

r(y, u, - U ) = ^ r ((y/o^) ( ^ +2)/2 , _ ^ + 2) , ^-(H-a)^ , (41) 



which is easily obtained either by directly inverting (p9|), or by properly changing the 
integration variable in (|i~3|). 



5.2 Parameter Duality Map 



The correspondence of the indeces involved in (39) 



2u 

M v : v — 



Z/ + 2 1 



is a one-to-one map of the semiaxis v £ (—2, oo) into itself with the only fixed point 
v = 0, so that (0, oo) is mapped into (—2, 0), and vice versa. Therefore for the case 
of the open universe it is enough to study cosmological solutions with either positive or 
negative v, getting then the results for the opposite sign of this parameter immediately 
from (|39|) and (|40|); we are intensively using this approach below. 

Moreover, it is easy to see that M. v is a duality map, that is, M v = M. v ~ l , or 
Mj 1 = Z, X being the unit map. The whole range of values of v is thus represented as 
a set of pairs of values reciprocal under the map A4 U , with v = reciprocal to itself. In 
particular, one reciprocal pair is v = 2 (relativistic gas) and v = — 1; it is straightforward 
to see that A 2 (kio) from @ and A^ 1 (kuj) from (||) satisfy (|D|). 

Note that for the physical parameter w entering the equation of state, the corre- 
sponding duality relation is 

w + 5/9 

M w : w — ► w>-l. 

w + 1 

It gives a one-to-one map of the interval (—1, —1/3) into the semiaxis (—1/3, oo) and 
vice versa; its only fixed point w = —1/3 corresponds to the Einstein quintessence, 
which once again separates one group of cases from the other. The reason for that is 
well known, and it is seen from the equation for the universal acceleration, which is 
obtained by the obvious manipulation with the governing equations (^): 

-31? = ~ 47T (P + 3p) + A = -Att(3w + l)p + A . (42) 

Therefore matter with w < —1/3 enhances the accelerating effect of vacuum, while 
matter with w > —1/3 diminishes it. It is only the Einstein quintessence, which is thus 
said to have no gravitational mass, that does not contribute to the acceleration. So, 
given the matter abundance M and k = —1, an expansion with w > —1/3 and certain 
A can in principle be related only to some expansion with w < —1/3 and smaller A, 
as it occurs according to (|39|). Of course, the Einstein quintessence solution cannot 
reciprocate with any other having a different equation of state. 
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5.3 Open Universe with u; ^> 1 



As a first example of using the transformation property, let us consider the asymptotic 
behavior of the expansion amplitude for the open universe (k = —1) when uj is large, 
which is easily found from fllpp and (p0[). Indeed, in this case uj~( v+2 ^ 2 — ► 0, and we 
have: 

A v {-u) =u 1 ' 2 

This result agrees with all the exact solutions (^2|), (I24T), and (^6|), as well as with the 
estimate (^9|), 

Thus the main term giving the large-time behavior of the open universe with uj 3> 1 
does not depend on the equation of state. The corresponding full asymptotic series, which 
proves to be converging, is found in section 6.5; all higher order corrections involved in 
it do depend on u, i. e., on the state equation. 



A -2v/(v+2) (0) 



1 + ^(1)] 



[1 + 0(1)] 



(43) 



6 Expansion Amplitude at cu Near the Cutoff Value 

As explained in section 2, cosmological expansion remains infinite, and its growth am- 
plitude is thus well defined, up to w = uj v for v > and closed universe (k = 1), which 
value forms thus a physical cutoff, in this case. For —2 < v < and open universe the 
solution exists for all positive values of kuj, so kuj = should be also treated as a cutoff. 
The expansion amplitude behavior near the cutoff value, which cannot be smooth, is of 
interest in both cases. 

We first consider v > 0, k = 1, and uj near uj u , that is, uj v — uj — > +0. In this case the 
main contribution to the integral ( |l5| ) for <j) u {ku)) comes from the vicinity of the point 
x = y u , where 

x~ v + x 2 - kuj ai u v - uj + (v + 2)(x - y v ) 2 . (44) 
Using this approximation it is not difficult to find that 

(j) u (u) = 1 ln(av - uj) + O(l) , k = 1, uj — >■ uj v - , 

+ 2 

and 

A„(w) = A u (0)K u (uj u - uj) 1/V ^ [1 + o(l)} , fc = l, (45) 

Calculation of the value of the constant K v requires more tedious analysis carried out 
in section 6.4, which provides 



Ln=l 



(46) 



where the coefficients C n {v) = OiAjn ), n — > oo, are given explicitly in (p4[). 

For Einstein quintessence, v = 0, and any curvature, approximation ( [i"4"|) works 
uniformly in > only for the function and its first derivative, but not for the second 
one. Therefore, with v formally set to zero, it gives a wrong coefficient at the quadratic 
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term (2 instead of 1). So for v = the right answer is given not by ([44|) with v = 0, but 
by the exact formula ( p4|) (rederived, by the way, in a different manner in section 7.3). 

Finally, when —2 < u < 0, we should consider the open universe only, k = —1, in 
which case we can write, for u — > +0: 

1 1 



4>v(kuj) 



Vx~ u + x 2 +uj \f. 



-l/2v 



1 



X u + X 2 
1 



+ ^-(2+^/^2 + X y^-z, + ^-(2+z/)/^2 

[recall that -l/2i/, -(2 + > 0]. So, by pi) and (Efl), 



-l/2i/\ 



M-o>) = [i + o(ar 1/2v )l , 4,(0) = 2~ 2 ^ +2 ) 



+0, 



(47) 



The exact result (|26| ) for yl_i(£;u;) agrees with this. Note that the remainder term in ([47|) 
diverges at to — > +0, hence so does the derivative of A J/ (— w), as expected. Thus, in fact, 



the behavior of the growth amplitude at both cutoffs, uj 



v > 0, and u = 0, v < 0, 



is similar, with the only difference that the limit value of the amplitude is zero in the 
first case and positive in the second. 



7 Power series in ku for A v (ku) and its implications 

7.1 Power series in ku for v > 

Repeated differentiation of the equality ([u]) in kuj provides 

[ln^(M] (n) = -4 n) (M, n = l,2,...; ^ . (48) 

Also, repeated differentiation of ( |T7| ) shows that, for i/ > 0, 

&\0)=(l) f°- ^172, n = l,2,.., (49) 

where the Pochhammer symbol (a) n is defined in a usual way, (a)o = 1, (a)n = 
r(a + n)/T(a) = (a + n — l)(a + n — 2) . . . (a + l)a, n = 1,2, and T(a) is the 
Euler gamma-function. [Note that the integral in (^) diverges at the lower limit when 
v < 0]. Therefore 

^U-f^tM", (50) 

n=l 



In 



4,(0) 

and the coefficients of this series are found to be (see Appendix) 

^ .«M, „. 1 , 2 ,.... (51) 

n! 2^/7r n 

The radius of convergence and other properties of this series are studied below. Since the 
radius proves to be non-zero, function A v {kiS) is regular at kcu = 0. That means that, 
for a given positive u, the expansion amplitude for the open universe is an analytical 
continuation of that for the closed one, and vice versa. 
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7.2 Case u = 2 

The result (|22| ) for v = 2 is easily obtained from the series @. Indeed, by © 

we have 



4 W) (Q) 

n! 



! r(§ + i)r(f + i) _ ! ^Fr(n + i) 1 



n 



2n+l ' 



nr(n + l) 20rn2™r(n + l) 

where we have used the double argument formula for the gamma- fucntion, 2 2 *- 1 r(z)r (z + 1/2) 
0FT(2z), with z = (n + l)/2. Thus 



In 



Ao(M 



4>(0) 



1 ~ (fc W /2) r 



n=l 



77. 



i ln(l - fcw/2) , 



and, according to (^0|), it is exactly the formula (p4|): A2(ku>) = 0.5(2 — /cu;) 1 / 2 . Note 
that the radius of convergence of the series ( |50"| ) coincides with the cutoff value for this 
case, 2 = W2- 



7.3 Case v = 

In this case formula ( ]5l| ) gives immediately </>i n ^(0) = l/2n, so 



In 



A (ku)) 



L MO) 



1 ^ (kuj 



n=l 



77 



2 ln (! - M ) 



that is, by (|2C|), AoikuS) = 0.5(1 — fccj) 1 / 2 , which is exactly the expression (^ij) . Here 
again the radius of convergence of series ( |50"| ) coincides with the cutoff, 1 = u;o, hinting 
that this should probably be a generic property. 



7.4 General case /y > 

The general result for v > comes from studying the behavior of the coefficients of the 
series (|50|) with large numbers n. As demonstrated in the Appendix, 



4 n) (0) 



1 1 



71! 



1 



l + O 



n — > oo . 



(52) 



+ 2 n \uj u , 

This means that the convergence radius of the series (|5^) always coincides, as one would 
now expect, with the cutoff value uj u , and allows us to rewrite the series as 



In 



A u (0) 



1 



v^T2 



71=1 



7? 



n=l 



1 



ln(l - fcw/w,,) + >T C n (i/) {ku/u) v ) n , 



(53) 



where, in view of (ml) and (52), 



1 1 



vV + 2 n 



2 



+ 2 



7T 



2ii 



n=l 



Oh, 



77- 



n — > oo . 



(54) 
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Expression ( p3[ ) is useful for computing A u (ku) in the whole range koj < uj v . Also, for 
k = 1 it determines completely the asymtotics of the amplitude near the cutoff lo u , whose 
main term is given above in (pj), (p3|). 

7.5 Series in inverse powers of u for — 2 < z/ < 0, k = — 1 

To obatin some series for A u (kuj) with negative ^ (and thus A; = —1), we combine the 
above result (|5^) with the amplitude transformation property (|40|) . Writing for brevity 
/i = —1vj{y + 2), the latter provides 



1 ,,1/2 1 

lnA y (-(j) =lno; 1 / 2 +^^lnyl At (-cJ^ /(iy+2) ) = In — - + . In 
VZ/ + 2 2 yzv + 2 



Mo) 

where we used the value of ^4^(0) from (|20|). Since, by definition, fi > 0, we can now 
replace the last term by its series expansion ( |53[ ) provided that uj~ v K u+2 ) < u;^, or 

uj > (55) 



[recall (y + 2), (— i/) > 0]. Under this condition we thus obtain the desired series for 
— 2 < < 0, k = — 1 in the form [the series coefficients are given in (1541)]: 



£.1/2 ! oo 

lnA„(-w) =ln— - + -== In [1 + «„(««;)] + y" C n (/i) M")]" , (56) 
2 V" + 2 n=1 

2^ . , -v/2 
U = , < k v (uj) = -. oWo < 1 . 



This series converges if condition ( |55| ) holds, and also provides the asymptotic formula 

,1/2 r 



A,(-u/) 



or 



*\ 2 



l + ( U <t>+W 



u -> oo, -2 < z/ < 0, fc = -1 , (57) 



which sharpens expression ( |43| ) by giving the exact order of the remainder. 

It is worthwhile to note that similar power series representations can be obtained 
also for the cosmological solution (|l3|), r(y, u, kui). However, the y-dependent coeffi- 
cients of these series are no longer expressed via T-functions, but, instead, through the 
hypergeometric functions. 



8 Discussion 

8.1 Cosmological Applications 

We have determined various significant properties of the cosmological solution, particu- 
larly, those of the expansion amplitude, and provided different expressions for the latter. 
The plots of A v {kui) as a function of kui for several 'popular' values of the parameter z/, 
including positive, zero, and negative ones, are given in Fig. 2. 
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An interesting question is what information on the universe can be gotten if the 
value of the expansion amplitude, say, A*, is known? Even if no other information is 



available, but A* > 1, from ( p0| ) and the monotonicity of the amplitude established in 
sec. 2, one immediately concludes that the universe is open, since A v {ku) < 1 for k = 1 
(and, of course, v > 0). If, on the other hand, < A* < 1, but the equation of state, i. 
e., parameter v, is more or less known, then equation 

A v {kuS) = A* 

provides, again due to the monotonic property, a unique solution for kto, which gives the 
sign of the curvature, and a relation between M and A [see (§)]. 

One way of inferring the amplitude from observational data is as follows. Formula 
d|]) for the density combined with (|7|) reduces to 8irp = A/y u+2 , therefore 



8tt P 



y 



u+2 



l/(w+2) 



(58) 



so that 

This equality holds at all moments of time. If, however, the expansion is in the asymp 
totic regime, ( |58|) and ( |l2|) allow for 



A v {kuj) 



l/(u+2) 



exp(- 



y(r) > 1 



(59) 



So, if the matter equation of state, the ratio Oa/Om, and the properly scaled [see (0)] 
cosmic age are known, the expansion amplitude is known as well, if the asymptotic stage 
of the expansion is reached. By (|58|), the latter requires 







A I 



l/{u+2) 



> 1 



(60) 



which, for v not much smaller than negative unity, means just Qa/^m ^ 1- But in 
case the dominant form of matter in the universe is quitessence with the equation of 
state close to the vacuum one, condition ( |60[ ) might hold, and the expansion be in the 
asymptotic regime, despite the fact that Oa/^m > 1 is still of the order of unity. 

For any moment of the expansion the estimates for the age of the universe come 
from combining expression (^) for the scale factor via observables with 'geometrical' 
inequalities (|33|): 



t < T F i(Q A /n M , v), k = -1; r > T F i(0, A /n M , v), k = l 



Here 



T F i(Q\/n M , v) 



u + 2 



hi 



1/2 



+ 1 + 



Oa 



M 



1/2 



(61) 



(62) 



is the exact age of the flat universe with the same parameters. If v, VL\/£Im, an d the 
cosmic age are known accurately enough, its comparison with the bound (^) allows for 
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an immediate determination of the curvature of the universe. If, instead, only the last 
two parameters and the sign of the curvature are known, the state equation parameter 
can be estimated. 

We now apply the above results to the currently favored cosmological model. It 
includes baryonic and dark matter, both with the same equation of state w = and 
total VLm = 0.3 ± 0.1, and vacuum with = 0.7 ± 0.1. The third component is 
relativistic particles (w = 1/3), but its abundance 0# < 0.001 is so small compared to 
the first two that it should be neglected. This results in exactly the studied model with 
v = 1 and Oa/^m = 1-5 -j- 4. Thus 

rik^.^r.!.^ ( 63 ) 



so, by (|58|), we are still rather far from the large time regime. 



According to (62) and (p3j), for our universe 

TFl = 0.7 -=- 1.0 . (64) 
The dimensional time is related to r by [see (R)] 



T = T^j-^ = T (0.7 ~r~ 0.8) x 10 18 sec = r (22.2 25.4) Gyr , (65) 
with 

A = S-npy = 8vr(0.7 ± 0.1)p c = (1.0 -f- 1.6) x 10~ 28 g/cm 3 

based on 

p c = (0.7 ± 0.08) x 10" 29 g/cm 3 , 
which corresponds to the range of the Hubble constant H = 70 ± 10 km/ sec ■ Mpc. By 



virtue of (64) and (pq), for our universe 

T Fl = tfi (22.2 -=- 25.4) Gyr = (15.5 + 25.4) Gyr , (66) 

that is, its age in case the universe is flat. Globular cluster data independently give the 
lower bound of the cosmic age as 

T u > (12 + 16) Gyr . (67) 

According to these data, our universe is open if T u = (12 -r- 15.5 — 0) Gyr, it could 
be either open or closed if T u = (15.5 div 25.4) Gyr, and only closed if T u > 25 A) Gyr, 
which is rather improbable. Therefore one can conclude that our universe most probably 
is open. 

Under this assumption we can now obtain a lower bound of the cosmic age, to check 
the consistency of our estimates. Indeed, using (^), we obtain an estimate of the 
parameter u> for our universe as 

(n A /n M f 3 

0J < ——o = m < 2 - 7 = ^est ■ (68) 

sum t u sinn t u 
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This and the basic expression (13) for the solution immediately produce 



r > 



(!iA/SJi;) 1/3 



dx 



> 



1.1 



dx 



or, by (H), 



y/x- 1 + X 2 +0J est Jo Vx' 1 + x 2 + 2.7 



T u > (10.0 -=- 11.4) Gyr 



0.45 , 



in a nice agreement with (|67[). 

Finally, to check the consistency of the current cosmological data, we can convert the 
estimate (61) for an open universe into an upper bound for the state equation parameter 
v, namely: 



v < 2 i T7 1 In 



V J 



1 



(69) 



By (|67D taken as the real range of the cosmic age, (|6q) and (p3|) , for our universe this 
provides 

v < 0.5 -=- 10 . 

This is almost entirely consistent with the chosen v = 1 (pressureless dust), except for 
a small range of parameters, when the ratio Q\ is close to its minimum observed value 
0.6 and simultaneously the cosmic age is close to its maximum of 16 Gyr. 



8.2 Some Generalizations 

In conclusion, we point out two possible generalizations of the present study. The first 
of them deals with the phase transition, that is, with a sudden change of the equation of 
state, when the value of w changes abruptly at some moment of time. This well 
as the one with a whole sequence of phase transitions, can be investigated in a fashion 
similar to the above. 

The second way to generalize the study is to consider multicomponent matter, when 

N N 

P=^2Pn, P=^2Pn, Pn = W n p n , W n > -1, U = 1, 2, . . . , N > 1 , 

n=l n=l 

and the components do not interact. That means that the conservation equation in (||) 
holds for every component separately, yielding p n = M n /a 3 ^ Wn+1 \ M n = const > 0, and 
thus the analog of the governing problem @ becomes 

N 

y 2 = y- Ul +y 2 -kco + J2 PnV' u \ 2/(0) = , 

ra=2 

with y, t and u normalized as in (0), (||) using w\ and M\, v n = 3w n + 1, and 
fi n > being the abundance of the nth component normalized appropriately. Writ- 
ing y(r, v\, v\, . . . , vn, ku>) for the solution, we have 

y(r, vi,vi,...,i/ N , kid) > y(r, uj, ku) 
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for any pertinent j; generally, adding every new component enhances the expansion. 
Assuming v\ is the smallest of all the powers, we see that at very large times the expan- 
sion is described by the one-component equation (||) with v = v\\ however, the growth 
amplitude depends essentially on the whole expansion history, in other words, on all the 
equations of state. 

As for the other results, all the estimates of sec. 4 remain true, including inequalities 
(|33|) for the ages of open, flat, and closed universe (of course, the expression for tfi 
depends now on all the parameters involved). The properly generalized transformation 



property (39) for the open universe solutions also holds; in fact, the number of such 



independent transformations is equal to the number of matter components, N. 
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Appendix. Analysis of Power Series of A v {ku) 

First we derive expression (|50|) for the derivatives ^^(O). By (|49|) , for n = 1,2,... we 
have: 

i/(2n + l) + 2 , 

i(n) fn\ r° X K2n+1)/2 dx 1 foo s^+T) 1 ds 



(0) 



r(n + l/2) ri/ w Jo (i + ^+2)»+i/2 u + 2Jo (l + s ) n+1/2 

The last integral is one of the known representations of the beta-function B(x,y) with 
the arguments 

u(2n + 1) + 2 vn 1 1 2n 

x = — — — — = + -, y = n + 



2(u + 2) u + 2 2 2 u + 2 

which is easily verified by means of the substituion s = tj {! — £). Therefore 

r(n + i) ^ + 2 + 2 2'i, + 2; 2nu + 2 r (n + |) 

i r(^ + i)r(^ + i 
2n r (n + \ 



and this is, in fact, expression (|50|). 

We now prove the asymtotic formula ([52]) for the coefficients of the seires ([5(]) in the 
general case. To do that, we use the large argument asymptoitcs of the gamma-function, 

T(z) = V2^ exp [(z - 1/2) In z - z] [1 + 0(1/ z)] , z -> +oo, 
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and convenient notations 



u + v = 1 . 



We thus write: 



z/ + 2' v + 2' 

4t\u) _ i r(«n + i)r( w + |) 



n! 2^n 

and then 

r (im + 1) r ( 



2ir exp(Q„) 



T(n + 1) 



.1 + O ( - 

n 



n — > +oo , 



(Al) 



(A.2) 



r(n + i) 

where 

Q n = (nn+1/2) ln(un+l)-(un+l) + (W) ln(tm+l)-(vn+l/2)-(n+l/2) ln(n+l) + (n+l) 

(un + 1 /2) ln(im) + (to) ln(wre) — (n + 1 /2) In n + 0(1 /n) = 
n ]n(u u v v ) + (1/2) lnu + O (1/n) . 
Intorducing this into (A.2), and the resulting equality into (A.l), we end up with 



1 



^ n) (Q) _ 

n\ 2^n 



V2^L(u u v v ) n 



1 + 



i (u u v v y 



y/u + 2 n 
It remains only to notice that 

2 \ 2/(^+2) 



i + o( - 

n 



71 — > +00 



(A3) 



+ 2 



z/ + 2 



i/ + 2 uj u 



in order to see that (A. 3) is identical to (52). 
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Figure 1: Parameters for the Closed Universe Solution 
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Figure 2. A u (km) for nonnegative v. 
- v = 0, 1 - v = 1, 2 - i/ = 2, AS - lower bound 0.5 A /(-A;a;) 
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Figure 3. A u (km) for negative u: 
\-v = -1/2, 2 - v = -1, 3 - v = -3/2, AS - lower bound 0.5 A /(-fcw) 
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